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rn ' Abstract 
vO ' 

The threshold transverse response functions Rx{q-,oj) for 3 He and 3 H are calculated using the 

O 

\q . AV18 nucleon-nucleon potential, the UrbanalX three-body force, and the Coulomb potential. Final 

o, 

Q^ states are completely taken into account via the Lorentz integral transform technique. Consistent 

o 

two-body 7r- and p-meson exchange currents as deduced using the Arenhovel-Schwamb technique 
^ ■ are included. The convergence of the method is shown and a comparison of the corresponding 

h : 

MEC contribution is made to that of a consistent MEC for the meson theoretical r-space BonnA 
potential. The response Rt is calculated in the threshold region at q=174, 324, and 487 MeV/c and 
compared with available data. The strong MEC contributions in the threshold region are nicely 
confirmed by the data at g=324 and 487 MeV/c although some differences between theoretical and 
experimental results remain. A comparison is also made with other calculations, where the same 
theoretical input is used. The agreement is generally rather good, but leaves also some space for 
further improvement. 

PACS numbers: 
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I. INTRODUCTION 



In a recent publication [l| the LIT technique 



was used to compute the transverse 



response Rt{o,uj) for 3 He with a quasi-modern potential. In that work we employed the 



configuration space BonnA potential [5j together with the Tucson-Melbourne' (TM') [(| NNN 
potential and the Coulomb force. Rather detailed numerical checks on our computational 
methods have been performed. As stated in [l( the BonnA potential was chosen because, 
being of boson-exchange character, meson exchange currents could be determined relatively 
uniquely. Although we had computed Rt{q-, ^) for a wide range of momentum transfers q and 
excitation energies u our choice of the BonnA potential did not allow a detailed comparison 



jj, |8|. That is because Golak et al [7| used the Argonne 



with other recent theoretical work 

n 

VI 8 (AVI 8) H NN potential with the UrbanalX (UIX) NNN potential [l0( while Deltuva 
et al [8j employed a coupled channel CD-Bonn [rJJ+A calculation. Despite these differences 
there is a great similarity between our BonnA results and those of 0, Q] ■ In particular the 
effects of meson exchange currents appear very prominently in the near threshold region of 
the response function. Hence the purpose of the present paper is to examine the threshold 
region again but this time with the AV18+UIX combination and with the addition of the 
calculation of Rt(q,uj) for 3 H as well. These results should be directly comparable to the 
results of |7j] except for the fact that with the LIT method we are able to consistently include 
the Coulomb interaction in the initial and final states. 

As indicated above meson exchange currents are relatively straightforward to obtain when 
the NN potential is of boson exchange type. This is not the case with the AV18 potential and 
thus a prescription is required for the construction of consistent meson exchange currents. 
Several methods based on interpreting the isovector part of the potential as due to an 



effective it and p exchange have appeared in the literature 



12 



13 



14j . These methods are 



similar in principle and differ mainly in form. Here we choose the technique of Arenhovel 
and Schwamb [14J whereas Golak et al employed the method of Riska 121 ] . 



The next section provides a very brie 
tails can be found in our previous paper 
the method of 



review of the calculational technique. Full de- 

I 

1[. Following that we describe our application of 



14j to construct effective ir - and p-exchange currents for the AV18 poten- 



tial. Finally we present, discuss, and compare our results with experiment as well as other 
theoretical calculaltions. 
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II. THE METHOD OF CALCULATION 



The transverse response Rt which depends on the transverse nuclear current density 
operator J T is given by 

R T (q,u) = Y, Mo ^df(^o\4^^f)-^f\M^^o) S(E f -E + q 2 /(2M T )-co). (1) 

Here Mt is the mass of the target nucleus, \I/o and denote the ground and final states, 
respectively, while E and Ef are their eigenenergies, 

(/i-£ )*o = 0, (h-E f )* f = 0. (2) 

Also in Eq. (2) h denotes the intrinsic nuclear non-relativistic Hamiltonian which includes 
the kinetic energy terms, the 2N and 3N force terms. In the present work the 2N + 3N inter- 
actions are taken as the AV18+UIX+Coulomb potentials. By decomposing the transverse 
current into electric and magnetic multipoles the response function itself can be written as 
a sum of multipole components via 

4-7T 

Mfc") = 2jXi E E( 2J + 1 )(^ A ( 3 ) 



2 Jo + - , , 

A=el,mag Jj 



where 



(RtYj X = J£ df{ q %\^ } {J,M))^ f {J,M)\q% [ )5{E } -E Q - u ), (4) 
J and M are the final state angular momentum and its projection, and Wjm) is given by 

I^m) = [^ A ®I*o(Jo))]^. (5) 

In Eq. ([I]) M is arbitrary while in Eq. ([5]) Tj* are the standard electric (A=el) or magnetic 
(A=mag) multipole operators. Note that in [l| we used two forms of the transverse electric 
operator, one containing only the current operators as used here, and another called the 
Siegert form which has a piece depending on the charge density operator as well. There it 
was found that in the threshold region either form gave identical results. Therefore here 
we use non-Siegert form of the transverse electric multipole operator. Further we note that 
nucleon form factors are the same as those used in l|, except for the neutron electric form 
factor which is taken from [151. 
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The techniques we use in calculating the response have been largely set out in 

ii a. 

Briefly, the Lorentz transform of the partial response (Rt) j A is given by 

The sum in ([6]) corresponds to transitions to discrete levels with excitation energy u n . 
In our A=3 case there exists only one discrete contribution corresponding to Ml elastic 
scattering. In (jSJ) the response is supplied with an additional superscript a. It specifies 
separate contributions to the response (R T yf of Eq. (@D,e.g. a given a determines the isospin 
of the final state. In addition it specifies contributions that correspond to components of 
the multipole operators with different nucleon form factor dependencies (for further details 
see These transforms are determined dynamically from 

$j A '°(<?, <tr, 0j) = $j X m$jm), I^jm Q > = [h-<r R + ^iY 1 Wj X m)- ( 7 ) 
Once computed these &j X,a (q,<jR,<Ti) are inverted separately to obtain (R T yj X ' a (q,u) and 

47T 



then R^' a from 



R^u) = — - T E( 2J + l)(^T)S A '°(g^)- (8) 



III. ONE-BODY AND TWO-BODY CURRENTS 

The one-body current in the present work consists of the non-relativistic current used in 
[l| plus all the relativistic corrections up to order M~ 2 i.e. it includes the spin and convection 
current terms which are of order M" 1 plus all the terms of order M -3 . We have calculated 
this operator from the expression for the corresponding single-particle matrix element of 
the form (p/| J|p^) given in 

Below we give a short description of the method of Arenhovel and Schwamb which 
we then use to obtain tt- and p-exchange currents for the AVI 8 potential. One begins with 
the well-known forms for the static nucleon NN potentials due to single tt or p exchange, viz 

K(r) = V^(r 1 -T 2 )(a 1 -V)(a 2 -V)J m7V (r), (9) 
V p {r) = V p (n • r 2 )(<7i x V)(tr 2 x V) J mp (r) (10) 

where 

p—mr 

J m {r)= • (11) 



If there were a continuous set of 7r-like exchanges then J mn (r) would be replaced by 



V*J m „ -> J n = V n / dm g^m) J m+mn (r) 

'o 



(12) 



which can be interpreted as a superposition of other PS-mesons with mass m+m,,. and 7rNN 
coupling constant densities of V^g^m). Similarly for the vector bosons one writes 



poo 

J P = V* dm g p {m) J m+m ,(r) 
Jo 



(13) 



where again the hypothetical vector mesons of mass m+m^ have pNN coupling constant den- 
sities of V 7r 5'p(m). Note that the masses of the hypothetical vector bosons here extend from 
m^ to oo. Application of the derivatives in Eq. gives the central and tensor potentials as 



1 f°° 
V (r) = -V„Jm*(r) / dm[g w (m) + 2g p (m)](m + m 7T ] 
<J Jo 



2 ^—mr 



1. 



V T (r) = J m A r ) J dm [gA m ) - 9p{ m )\ ( m + m 7v) 2 F T [(m + m^r] e 



(14) 



where 



Fximr) 



1 + — I 1 

mr 



1 

mr 



and where a delta function 8(r) is removed from V (r) by imposing the condition 



dm[g w (m) + 2g p (m)) 



The change of variable 



m 



s tan 



7T 



(x + 1) 



and the introduction of an iV-point Gaussian integration in x=[-l,l] gives 



poo N 

/ f(m)dm -> wj f(m(xj)) 
Jo j=1 



where 



w j 



7T 2 

s— sec 
4 



Wj . 



(15) 



(16) 



(17) 



(18) 



Here the Xj and Wj are the N abscissae and weights for the quadrature integration and s 
is a parameter to be discussed shortly. If iV values of r are selected, say Tji = 1...N, then 
Eqs. (1141) can be solved for g n (mj) and g p (m,j) j = 1...N. We follow Ref. 14| in using their 
method to choose N values of r between r mir , = 0.01 fm and r mnT =12 fm . These values of 



the g n and g p thus give an exact fit to the AV18 values of V c (r) and V T (r) at the chosen 
iV values of r. The parameter s is adjusted to give the minimal absolute deviation between 
the original potential and its new parametrization when integrating the absolute deviation 
from r min to r max . Up to this stage the 5-function removing condition Eq. (IT51) has not yet 
been imposed. Ref. [14] imposes this condition by modifying the values of g,r(A) and g p {N) 
i.e. the contributions of the largest mass terms such that Eq. (Tl5l) is numerically satisfied. 
Thus if the least squares fitting yields parameter values such that 

/ dm[g n (m) + 2g p (m)} -> Wj^K) + 2g p {m j )} = C^O (19) 

J 3 

then the replacements 

9,{N) - g° ld {N) - and g p {N) = g° ld (N) - (20) 

Sw N F 6w N 

will satisfy Eq. ( fl5l) at the expense of spoiling the short range fit at r=0.01 fm. However 
the fit is still very good for r greater than approximately 0.06 fm depending on the choice 
of N. 

Finally the it- and p-meson exchange current multipoles obtained from the above method 
are, apart from slight modifications, equal to those listed in Appendix C of [1]. These 
modifications in the case of 7r-exchange currents are the following: (a) here we use 
/ 2 =0.075 and m,=0.70 fm" 1 , (b) H*(r) = £f=i w 3 g w ( mj ) J mn+m . (r), (c) $^(g,r) = 
YlJ=i'^'j9^( m j) 4 ><y a\{ ( li r i' m j)i where the functions 4>^j(q, r,m) are defined in [171 ] . The mul- 
tipoles of the p-exchange currrents are obtained from the 7r-meson exchange currents by the 
replacements: H n (r) — > H p (r), g w (rrij) — > g p {mj) and by inserting into each equation the 
factor 

fill, 

-6(-lW >■ (21) 

IV. RESULTS AND DISCUSSION 

For the calculation of R? we take into account electric and magnetic multipole transitions 
to final states with total angular momentum Jf up to 5/2 and in the case of final isospin 
Tf = 1/2 we include transitions to the J/=7/2 state in addition. The latter lead to a small 
contribution in the peak region for g=174 MeV/c (for this rather small q the " quasielastic" 
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peak is a part of the threshold region). The three-nucleon ground states and the various LITs 
are calculated using expansions in correlated hyperspherical harmonics [4]. The expansions 
in hyperradial and hyperspherical basis functions is made such that the convergence errors 
of the various LITs are considerably smaller than 0.5% in the energy region of main interest 
(-10 MeV < or < 20 MeV). As the resolution parameter 07 we choose 10 MeV, but for 
the two dominant threshold contributions (Ml and M2 transitions to Tf =1/2 states with 
Jf = 1/2 and 5/2, respectively) we take 07 = 2.5 MeV extending the HH expansion further 
in order to guarantee also in these cases the requested small convergence error. Before 
carrying out the inversion of the LIT we first subtract the elastic Ml contribution. The 
inversion itself proceeds in the following way. The above mentioned Ml and M2 transitions 
are inverted separately. Because of their different break-up thresholds the LITs of the 
remaining multipole contributions are individually summed up and inverted for the Tf = 1/2 
and 3/2 states. As to the inversion method we use our standard expansion of tie LIT oyer 



given basis functions, where also the threshold behavior is incorporated (see [J, [18|, Il9|). 
Note that in case of 3 He we also take into account the effect of the Coulomb barrier on the 
threshold behavior. 

In the first part of the discussion we consider some aspects of the consistent MEC for 
the AV18 potential. In Fig. 1 we show for 3 He the LIT of the MEC contribution for the Ml 
and El transitions to the Jf = 1/2, Tf = 1/2 final state at g=174 MeV/c for three values of 
the Gaussian integration parameter N discussed in section III. One sees that N = 12 does 
not yet lead to a convergent result, but that the N = 14 and N = 16 results are almost 
identical. The convergence of other analyzed MEC transitions looks very similar. Thus, in 
order to be on the safe side we use N = 16 for the calculation of the MEC contribution of 

Rrp- 

Next we make a comparison of the consistent MEC between BonnA and AV18 potentials 
(as mentioned before the unique BonnA-MEC was calculated by us in [l|). The MEC effect 
on the threshold response is much stronger for the magnetic than for the electric transition. 
Therefore we investigate the magnetic MEC contribution considering the LIT of the total 
MEC magnetic transition strength (sum of LITs of all MEC magnetic multipole transitions). 
In Fig. 2 we show the 3 He results for the two Tf channels separately at g=324 MeV/c. One 
sees that the shape of the LITs for BonnA and AV18 MEC strength is very similar, in 
particular for the Tf = 1/2 channel. In fact the main difference consists in the overall 
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strength. In comparison to the BonnA case the AV18 MEC strength is reduced by about 
12% (Tf = 1/2) and 9% (Tf = 3/2). For the other two q values, not shown in the figure, 
the shapes are also very similar, but here one finds somewhat different reductions of the 
overall strength, namely 9% and 15% (T f = 1/2) and 4% and 13% (T f = 3/2) at q = 174 
and 487 MeV/c, respectively. These results show that the ratio of the BonnA to AV18 MEC 
strength grows with increasing momentum transfer. 

Now we turn our attention to the comparison of our Rt results to the data of Retzlaff 
et al 20] (see Fig. 3). In case of 3 He one has a rather good agreement of theoretical and 



experimental transverse response functions for the two higher q values. The MEC contribu- 
tion is essential for reaching this agreement. At q=174 MeV, however, the theoretical Rt 
underestimates the data below 11 MeV. In the triton case the situation looks worse. Already 
for the two higher q values one finds a slight underestimation of the data, in addition the 
discrepancy becomes even larger at the lowest q. The relativistic contributions originating 
from the use of a relativistic one-body current are not negligible already at g=323 MeV/c 
and make the discrepancy theory-experiment somewhat larger. One can conclude that the 
present agreement between theory and experiment is not bad, but certainly not very good. 
It seems that a different nuclear force does not improve the situation, since our 3 He re- 
sults at g=174 MeV/c with the BonnA+TM' potential from [1] is almost identical to the 
AV18+UIX result (the 3 H case was not considered in [1]). Additional currents involving the 
A resonance, up to now only partially considered in the literature for the threshold kine- 
matics (see [a, l2l|), could probably lead to a small improvement. On the other hand, as to 
the experimental cross sections below the break-up threshold, one cannot exclude that the 
data are systematically a little bit too high. Thus, presently, one cannot speak of a serious 
disagreement of theoretical and experimental results. 

Finally we want to make a brief comparison of our results for the threshold response with 
other calculations with the same theoretical input. As already mentioned above there is 

Q 

the calculation of [7J, where a nonrelativistic one-body current and a consistent AV18-MEC 
have been taken as the current operator. For the small energy range up to the three-body 
break-up threshold there exists in case of 3 He another calculation [2l| . Here we only consider 
the results with a nonrelativistic one-body current, since the MEC of 21| includes also an 
additional current involving the A resonance, which makes a consistent comparison more 
difficult. In Figs. 4 and 5 we illustrate the various results (note no results available from [?| 
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with one-body current only). For the Rt of 3 H one has quite a good agreement between the 
different theoretical calculations. Some differences are visible below the three-body break- 
up threshold and above 14 MeV, only for the highest q value one finds also some small 
differences for other energies. For the latter a not unimportant part of the difference can be 
explained by the use of the three-momentum transfer squared for the nucleon form factors 
in case of 7|, while we take the four-momentum transfer squared. This leads to a relative 
reduction of the Rt of by a little bit more than 1% (note that at g=487 MeV one has a 
nuclear recoil energy of about 50 MeV). For the Rt of 3 He we first consider the result of 21] 
(see Fig. 5). Up to about 1 MeV above threshold our results are about 20% smaller. Close 
to the three-body break-up threshold differences are reduced to 9% (g=174 MeV/c), 3% 
(g=323 MeV/c) and 0% (g=487 MeV/c). Here we should mention that for the longitudinal 
response function Rl a much better agreement between the two calculations was obtained 
[31. For the 3 He results in Fig. 4 one has to take into account that there is no Coulomb force 

n 

in the final state interaction of |7(. This should explain why differences with our results are 
considerably larger than for the 3 H case, even though the effects seem to be a bit too large 
in order to be caused by the Coulomb force alone. 

Summing up we can say the following. For the AV18+UIX nuclear interaction we have 
calculated the (e, e') transverse response function Rr{q, oj) of 3 H and 3 He close to the break- 
up threshold at three momentum transfers (g=174, 323, and 487 MeV/c). Besides the 
one-body current a consistent isovector MEC for the AV18 potential has been employed. 
The consistent MEC has been constructed using the Arenhovel-Schwamb method, where 
an expansion over fictitious meson masses is made. The convergence of the expansion has 
been shown for our results. In comparison to the MEC contribution of the meson theoretical 
BonnA potential one obtains very similar results for the energy dependence, but the absolute 
size differs somewhat and is function of q and the final isospin channel. Relativistic effects 
for the one-body current have been included and lead to non negligible contributions for the 
two higher q values. A comparison with experimental data has been made and, as already 
known before in the literature, MEC contributions lead to a much improved agreement with 
data. As discussed in the text the comparison theory-experiment is not bad, but certainly 
not yet completely satisfying. The agreement with other theoretical calculations with the 
same theoretical input is overall rather good, but some differences, in particular below the 
three-body break-up threshold, are also evident. 
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FIG. 1: Lorentz integral transforms of the Ml MEC transition strength to the final state J n = l + /2, 
T f =1/2 (left panel) and of the El MEC transition strength to the final state J 7r = l _ /2, T f =l/2 
(right panel) calculated with a different number TV of hypothetical meson masses (see text): iV=10 
(dashed), N=12 (dash-dotted), iV=14 (solid), N=16 (dotted). 
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FIG. 2: Lorentz integral transform of the total inelastic magnetic MEC contribution for final 
isospin channel Ty=l/2 (left panel) and Tf=3/2 (right panel): consistent tt and p MEC for AV18 
potential evaluated for 3 He with AV18+UIX potentials (dashed), consistent ir and p MEC for 
BonnA potential evaluated for 3 He with BonnA+TM' potentials (solid), latter result renormalized 
with factors 0.88 (T/=l/2) and 0.91 (2)=3/2) (dotted). 
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FIG. 3: Comparison of experimental and theoretical results for the threshold 3 He (upper panels) 
and 3 H (lower panels) transverse response function Rt as function of internal excitation energy 
Wint at three momentum transfers q: 174 MeV/c (left panels), 324 MeV/c (middle panels), 487 
MeV/c (right panels). Theoretical results with different current operators: relativistic one-body 
current (dashed), relativistic one-body current + MEC (solid), nonrelativistic one-body current + 



MEC (dotted). Experimental data from 



20] 
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FIG. 5: Comparison of Rt with nonrelativistic one-body current from present work (solid) and 
from [2l| (dashed) at various q as indicated in figure (note that the result at q=2A7 fm _1 is 
multiplied by a factor of 5). 
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